中国空间探测领域40年发展

中国空间科学学会成立的40年，是中国空间探测逐渐走进世界舞台的40年，空间探测极大推动了空间科学和相邻学科的发展，也影响到经济、军事和日常生活诸多方面.本文简要回顾了从空间探测专业委员会成立的1980年至今，中国空间探测领域的主要发展历程，包括探空火箭、高空气球、科学卫星、月球与行星探测、载人航天空间探测、遥感卫星地面站等主要项目、进展和所取得的成果，对未来若干年空间探测的发展进行了展望.      

Assessment of the close approach frequency and collision probability for satellites in different configurations of large constellations

The Earth orbital environment is drastically changing due to an intensification of the space activities. In particular, several projects of large constellations, proposed for the next years for communications purpose like global internet access, Internet of Things, or for Earth observations, will lead to the deployment of several thousands of new satellites at an unprecedented rate. It is a crucial challenge for space traffic management, which will deal with a great number of satellite conjunctions, potentially causing a collision with damaging consequences for the constellation itself and the space environment sustainability.In this paper, we investigate the close approach frequency and the cumulative collision probability for each referenced constellation. For this purpose, we compute the orbital evolution of satellites in different constellations during the lifecycle, from the deployment to the decommissioning, and we apply the CUBE algorithm and the Foster method to assess the collision probability with the background space debris population assuming a constant uncertainty in position. We show the variation of risk defined by the close approach frequency and the cumulative collision probability as a function of the proposed configuration. In particular, satellites of the Iridium and Kuiper constellation, but also satellite of the Telesat constellation on polar orbits are the most exposed at a collision. Moreover, the decommissioning phase contribute for a major part to the final cumulative collision probability.     

A comprehensive survey on the methods of angle of attack measurement and estimation in UAVs

Angle of Attack (AOA) is a crucial parameter which directly affects the aerodynamic forces of an aircraft. The measurement of AOA is required to ensure a safe flight within its designed flight envelop. This paper intends to summarise a comprehensive survey on the measurement techniques and estimation methods for AOA, specifically in Unmanned Aerial Vehicle (UAV) applications. In the case of UAVs, weight constraint plays a major role as far as sensor suites are concerned. This results in selecting a suitable estimation method to extract AOA using the available data from the autopilot. The most feasible and widely employed AOA measurement technique is by using the Multi-Hole Probes (MHPs). The MHP measures the AOA regarding the pressure variations between the ports. Due to the importance of MHP in AOA measurement, the calibration methods for the MHP are also included in this paper. This paper discusses the AOA measurement using virtual AOA sensors, their importance and the operation.     

Influence of thickness and initial groove angle in M–K model on limit strain of 7B04 by considering through-thickness stress

Marciniack–Kuczinski (M–K) model is widely used to predict material’s forming limit curve (FLC). The prediction of FLC traditionally neglected through-thickness normal stress. However, it cannot be neglected in some forming processes. Much work has been done to study the effect of through-thickness normal stress on FLC with constant through-thickness normal stress or constant ratio of through-thickness normal stress and maximum principal stress. In addition, based on Nakazima test process, the ratio of through-thickness normal stress and maximum principal stress has been derived, which was a function of instantaneous thickness and loading path. Here, initial groove angle in M–K model was not considered. In this paper, uniaxial tension tests and Nakazima tests were performed on 7B04 aluminum alloy. Based on Hill 48 yield criterion and M–K model, the prediction model of FLC was established. The increase of thickness can enhance FLC. Meanwhile, it is necessary to consider through-thickness normal stress and initial groove angle in prediction model. On the left side of FLC, the effect of initial groove angle on FLC is weakened by increasing sheet thickness. On the right side of FLC, the effect of initial groove angle on FLC is strengthened by increasing sheet thickness. On the right side of FLC, the relation between limit strain points with different thicknesses is linear under one certain loading path. Thickness has decisive effect on through-thickness normal stress level and the changing trendy of through-thickness normal stress during calculation is different under different stress condition.     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

导航卫星中断频次的分析方法

面向导航卫星中断频次的定量分析需求,该文分析导航卫星中断产生的主要原因,给出中断频次分析的流程,并针对分析过程中的三个关键问题,研究提出具体实施方法,包括通过相关性分析快速定位底层中断事件,通过中断树建立指标分析模型,并融合在轨数据、地面试验数据快速预估得到底层功能异常率等。最后通过示例进一步说明中断频次分析过程。该文方法已应用于北斗导航卫星工程。     

Orbit-attitude-vibration coupled dynamics of tethered solar power satellite

A new orbit-attitude-vibration coupled dynamic model of the tethered solar power satellite (Tethered SPS) is established based on absolute nodal coordinate formulation. The Hamilton’s equation of the system is derived by introducing generalized momentum through Legendre transformation. The correctness of the proposed model is verified by an example. The dynamic characteristics of the Tethered SPS are studied using the symplectic Runge-Kutta method. Simulation results show that the orbital radius and the total energy of the system are well preserved. The attitude of the system is unstable when the mass of the bus system is small. However, the attitude stability is dependent on some other parameters of the system, which requires further studies. It is also found that the average tether force/deformation can be roughly estimated by simplifying the solar panel as a particle. The proposed model can be used to study the orbit-attitude-vibration coupled dynamics and control problems.     

先进战斗机过失速机动大气数据融合估计方法

可控的过失速机动是先进战斗机超机动性能的重要标志，飞机飞行包线的扩大已超出传统的大气数据系统测量范围，可靠的迎角、侧滑角、总压、静压等飞行大气数据是制约先进战斗机过失速机动中飞行控制的关键因素。以中国推力矢量验证机为对象，基于过失速机动飞行试验的数据，开展大气参数估计与验证研究。结合过失速机动的时间与空间特性，研究了基于风速、地速、空速矢量和惯性姿态、导航参数的大气参数融合计算方法；针对过失速大迎角状态下飞机周围气流非定常、模型非线性导致的融合大气参数误差的复杂特性，进一步构建深度神经网络，对机动状态融合迎角、侧滑角的强非线性误差进行拟合。仿真和飞行试验表明：该方法可在大迎角飞行状态下实现主要大气参数的融合估计，过失速机动过程中融合迎角误差优于2.3°，融合得到的大气参数可为过失速大迎角机动飞行控制提供可靠的大气参数状态反馈。     

叶片式预旋喷嘴出口流场实验研究

叶片式预旋喷嘴具有尺寸小,落后角大的特点。为了详细研究小尺寸预旋喷嘴的预旋性能,采用五孔探针对叶片式预旋喷嘴的出口流场进行了实验研究。测量了Ma=0.2,0.3时喷嘴出口的压力分布、速度分布和出口气流角度分布,实验获得了喷嘴的落后角和预旋效率,并进行了与实验工况相同的数值计算。通过实验获得的总压云图以及速度云图,可以发现叶片式预旋喷嘴的端壁二次流损失、尾迹损失严重,有明显的边界层分离现象。Ma=0.2时,喷嘴Re数为5.76×104,落后角2.84°,实验测得的预旋效率为0.73;Ma=0.3时,喷嘴Re数为1.06×105,预旋效率提高至0.77。实验模型端壁的影响使预旋效率实验结果偏低6.5%左右。数值结果与实验测得各参数符合较好:数值结果与测得的喷嘴出口截面平均总压、静压偏差在1%以内;出气速度、周向速度以及出气角度与实验结果偏差在4%以内。数值计算表明,叶片式预旋喷嘴的预旋效率基本不受压比影响,随Re数增大先增大后基本不变,最后基本稳定在0.85。     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.